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Abstract 

The spin-singlet P-wave state of charmonium, h^Pi), has been observed in the decay ip(2S) — > 
ir°h c followed by h c — ► 77? c . Inclusive and exclusive analyses of the M(h c ) spectrum have been 
performed. Two complementary inclusive analyses select either a range of energies for the photon 
emitted in h c — ► 7ry c or a range of values of M(i] c ). These analyses, consistent with one another 
within statistics, yield M(h c ) = [3524.9 ± 0.7 (stat) ± 0.4 (sys)] MeV/c 2 and a product of the 
branching ratios B f (ip{2S) -► vr°/i c ) x B fc (/i c -► 7%) = [3.5 ± 1.0 (stat) ± 0.7 (sys)] x 10~ 4 . When 
the r) c is reconstructed in seven exclusive decay modes, 17.5 ± 4.5 h c events are seen with an 
average mass M(h c ) = [3523.6 ± 0.9 (stat) ± 0.5 (sys)] MeV/c 2 , and B^B h = [5.3 ± 1.5 (stat) ± 
1.0 (sys)] x 10 -4 . If combined, the inclusive and exclusive data samples yield an overall mass 
M(h c ) = [3524.4 ± 0.6 (stat) ± 0.4 (sys)] MeV/c 2 and product of branching ratios B$B h = [4.0 ± 
0.8 (stat) ± 0.7 (sys)] x 10~ 4 . The h c mass implies a P-wave hyperfine splitting AMhf(I-P) = 
(M(1 3 P)) - M(l 1 Pi) = [1.0 ± 0.6 (stat) ± 0.4 (sys)] MeV/c 2 . 

PACS numbers: 14.40. Gx, 13.25.Gv, 13.20.Gd, 12.38.Qk 
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Since the discovery of the J/ip, the first bound state of a charmed quark c and charmed 
antiquark c 0,0], the cc (charmonium) spectrum has provided many insights about quarks 
and the forces holding them together. The charmed quark was the first to be found with a 
mass larger than the characteristic scale of quantum chromodynamics (QCD). Charmonium 
bound states thus could be treated starting from a nonrelativistic description 0]. One 
could calculate decay rates and level splittings and thereby determine the magnitude of the 
strong coupling constant as at the charm mass scale, and the Lorentz structure of the force 
confining quarks (see, e.g., (HIHEf for reviews.) 

The hyperfine (spin-spin) splittings in charmonium S-wave states are appreciable 0, H[ : 

AMhf(IS) = M(J/ij) - M(r} c ) ~ 115 MeV/c 2 , 

AM HF (2S) = M(^(2S)) - M(r]' c ) ~ 49 MeV/c 2 . (1) 

For an interquark potential V(r) = Vs(r) + Vy(r), the sum of vector Vy(r) and scalar 
Vs(r) contributions, only the vector part contributes to the spin-spin splitting [i. IH M. Eo|. 
giving rise in lowest order of l/m c (m c is the mass of the charmed quark) to a spin-spin 
interaction perturbation 

VSS{V) = Mr) = 9ml 6 (r) - (2) 

The second equality on the right-hand side is obtained when one takes Vy{r) = 4a 5 /(3r) 
and neglects the slow variation of as with scale. The resulting local spin-spin interaction 
then contributes only to splittings in S-wave states. Taking account of the scale dependence 
of as 9, 10] and Xcj wave function variations, one finds at most a few MeV/c 2 splitting 
between the \ l P\ state h c and the spin-weighted average (M(1 3 P)) of the 3 Pj states \ c j Q: 
(M(1 3 P)) = [M(l 3 P )±3M(l 3 Pi)±5M(i 3 P 2 )]/9 = (3525.4±0.1) MeV/c 2 . Small splittings 
AMhf(IP) = (M(1 3 P)) — M{l l Pi) are also consistent with a wide variety of estimates in 
potential models [ll[ and non-relativistic QCD 12], as well as with lattice gauge theory 
estimates [l3|. Values of |AM H f(1P)| larger than a few MeV/c 2 could indicate unexpected 
behavior of the vector potential Vy(r), unexpectedly large distortions of the masses of the 
l 3 Pj = Xcj states due to coupled- channel effects, or - in lattice theory - effects of light-quark 
degrees of freedom. 

The low-lying charmonium spectrum is illustrated in Fig. 1. The Xcj can be easily 
populated by radiative transitions from the il){2S). Their subsequent radiative decays to 
J/i/j also are prominent. In contrast, the h c = l x Pi cc state is not easily produced. It can 
be produced in the pp direct channel, and a few events were seen at the CERN Intersecting 
Storage Rings (ISR), clustered about M(h c ) = 3525.4 ± 0.8 MeV/c 2 Q. The significance 



of the signal was 2.3cr. Stronger evidence was presented by Fermilab Experiment E760 in 
the channel pp — >• h c — > 7[ J/ip 15], with a combined branching ratio 



x 10~ 7 (3) 



pp — > IL C — > 71 J l/J [lOIJ ■ W1LI1 A, CUIUUlIieU Ul clIlCIllIlg I clLlU 

(1.7 ± 0.4) x 10" 7 < B(h c -> pp)B(h c -> vr J/ijj) < (2.3 ± 0.6) 

for M(h c ) = 3526.2 ± 0.15 ± 0.2 MeV/c 2 , with an additional possible shift of up to ±0.4 
MeV/c 2 due to resonance-continuum interference. However, E835, the sequel to E760 with 
three times its integrated luminosity, did not confirm the E760 signal [lfj, E3] • Instead, a 
signal with ~ 3<x significance for pp — » h c — > •yr] c — > 777 was reported recently with 
M(h c ) = 3525.8 ± 0.2 ± 0.2 MeV/c 2 , width T < 1 MeV, and (10.0 ± 3.5) eV < T(h c -> 
pp)B(h c -»■ 77 c7 ) < (12.0 ±4.5) eV. 
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FIG. 1: The low-lying charmonium (cc) spectrum and some observed transitions. The bold-faced 
lines labeled "7r°" and "tei" denote the respective transitions ip(2S) — ► ir°h c and h c — > tei^c 
discussed in the present paper. 



The decay ip{2S) — > 7r /i c can occur via isospin mixing (e.g., 7r -^ mixing) in the neutral 
pion |l8[. Previous experimental upper limits on the branching ratio for this process are 
B^ = B(if>(2S) -> ir°h c ) < 42-80 xl0~ 4 for M(h c ) between 3500 and 3535 MeV/c 2 , and 
~ 15 x 10- 4 for M(h c ) ~ 3525 MeV/c 2 , where B fc = B(/i c -> 7r/ c ) El. Ko 
estimated B^ ~ 30 x 10 -4 . A recent theoretical range is £>,/, ~ (4-13) x 10~ 4 j2l|. 

The decay /i c — > 7r/ c is an electric dipole (El) transition whose matrix element should be 
the same as that for the decays Xcj —> jJ/i*- Estimates 2^1 of T{h c — > 7ei^ c ) range between 



160 and 560 keV; a recent value is 354 keV j23[. The hadronic and photon + hadronic decay 
rates of h c are not as well estimated, but the total width T(h c ) is generally found to be 1 
MeV or less, with Ref. [23] obtaining 0.94 MeV and hence B h = B(h c -> 7E1 r/ c ) = 37.7%. In 
other treatments this branching ratio can be larger; it is rarely smaller. In ip(2S) — > ir°h c 
the polarizations of the h c and ip(2S) should be almost identical, since the spinless tt° is 
expected to be emitted in an S wave. The subsequent El transition h c — > 777c should then 
lead to a photon with distribution ^(cos^) ~ 1 + cos 2 9 with respect to the beam axis. 

The present paper describes the identification of h c at the Cornell Electron Storage Ring 
(CESR), using the CLEO III and CLEO-c detectors, via the sequential process 

e + e~ ^^(25) (3686) ^ n°h c , h c - lElVc , vr -> 77, (4) 
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illustrated by the bold arrows in Fig. [T] labeled "7r°" and "tei," respectively. Exclusive 
reconstruction of r] c decays in seven modes permits observation of h c with convincing 
significance and little background, while inclusive analysis in which the r] c is not recon- 
structed provides a better measurement of M(h c ) and of the combined branching ratio for 
ip{2S) -> n°h c , h c -> jtj c . 

We mention relevant aspects of the CLEO detector in Section II. An overview of inclusive 
and exclusive analysis methods is presented in Section III. We then describe background 
sources and suppressions (Sec. IV), data sample and event selection (Sec. V), Monte Carlo 
samples (Sec. VI), the extraction of signal from the data (Sec. VII), and systematic errors 
(Sec. VIII). The combined results of the different analyses are presented in Sec. IX. A 
summary and discussion of the results are given in Sec. X. 

II. THE CLEO DETECTOR 

The data upon which the present report is based were taken with the CLEO III and 
CLEO-c detectors, described in detail elsewhere \24L l25l . l26l . |27| . Elements critical for the 
analyses presented here are the calorimeter and, for the exclusive analysis, the charged 
particle tracking and particle identification systems. The barrel (80% of Att) and endcap 
(additional 13% of 47r) electromagnetic calorimeters consist of a total of 7800 thallium-doped 
cesium iodide (Csl) crystals. Their excellent resolutions in position and energy (2.2% at 
E 1 — 1 GeV and 5% at 100 MeV) are a major source of sensitivity and discrimination against 
background in identifying the chain of decays ip{2S) — > 7i°h c — > 7r°7?7 c , and in measuring 
M(h c ). Pion/kaon separation is performed utilizing the energy loss in the drift chamber, 
dE/dx, and photons in the Ring-Imaging Cherenkov (RICH) counters. The combined dE/dx 
and RICH particle identification system has an efficiency of > 90% and misidentification 
rates of < 5% for both tt^ and K . Approximately one-half of the data sample used an 
upgraded configuration, denoted CLEO-c, with an inner drift chamber detector sensitive to 
longitudinal position [28j . 

III. OVERVIEW OF ANALYSES 

In the analyses described here one starts by looking for the neutral pion emitted in 
i)(2S) -> n°h c , expected to have an energy of E(n°) ~ 160 MeV for M(i[)(2S)) = 3686.111± 
0.025 ± 0.009 MeV/c 2 [29] when M(h c ) ~ 3525 MeV/c 2 , and the El photon emitted in 
h c — > IeiVc, with an expected energy in the h c rest frame of E^ei) — 502 MeV for M(t) c ) = 
2981.8 ± 2.0 MeV/c 2 |3(|. One takes advantage of the good energy resolution of the CLEO 
electromagnetic calorimeter by searching for an enhancement in the spectrum of masses 
M(h c ) recoiling against the 7r°, 

M(h c ) = [M 2 0j(2S)) - 2M(tfj(2S))E(iT ) + M 2 (tt°)] 1/2 , (5) 

reducing background by selecting a range of El photon energy E^tei) or r] c mass M(rj c ) in 
the transition h c — > 7?7 C , with 

M( Vc ) = {M 2 (h c ) -2E{ lE1 )[E(h c ) +p(vr )cos^(vr ,7 E1 )]} 1/2 . (6) 

Here E(h c ) and p(tt°) are the h c energy and the magnitude of the 7r° three-momentum in 
the ip{2S) rest frame, while #(7t°,7ei) is the angle between the 7r° and 7ei in that frame. 
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A search that is inclusive with respect to the i] c decay, i.e., one that imposes no further 
requirements on the r\ c decay products, exploits the full event yield. With a sample of 
approximately three million i/;(2S), an estimated product branching ratio B^Bh — 4 x 10™ 4 , 
and an estimated efficiency of about 15%, one expects about 180 counts in the h c peak in 
inclusive analyses, albeit on top of a background several times larger. 

An exclusive analysis, in which specific decay modes of the i] c are detected, benefits from 
much lower backgrounds with reduced efficiency. In the present analysis nearly 10% of all 
7] c decays are reconstructed, leading one to expect ~ 18 events with little background. The 
method is validated by reconstructing the more abundant 7] c decays in the direct reaction 
if>(2S) -> jr]c for which B(ip(2S) -> jr) c ) = (3.2 ±0.6 ±0.4) x 10" 3 BJ. (The Particle Data 
Group average of other measurements is (2.8 ± 0.6) x 10 -3 Q.) 

The following features are common to both inclusive and exclusive analyses. The sensi- 
tivity of the search for ip(2S) —>■ n°h c — > 7r°7?7 c depends upon the degree to which the 7r° 
peak can be recognized above a background which rises sharply as 7r° energy increases. Thus 
understanding of E(tt°) resolution is central to observation of the h c in this process. It is 
also crucial in pinning down the mass of h c . 

Because the signal ir° in ip{2S) ~^ ^°h c is expected to have fairly low momentum, its 
decay photons tend to be back-to-back in azimuth. Mismeasurements of their energies are 
partly compensated by the mass constraint used when combining them into a it candidate 
and thus affect the 7r° detection probability only minimally, resulting in a narrow distribution 
in 7T° energy and therefore in M(h c ), as will be seen in the specific analyses described below. 

At E 7 ~ 500 MeV (the energy of the expected signal for h c — > JeiVc), the experimental 
resolution of the photon energy is comparable to that expected from Doppler broadening of 
the h c when the photon is observed in the ip(2S) rest frame (~ 10 MeV). One can correct 
for this broadening using information on cos9(n°, 7ei) as in Eq. (JBJ). 

Two complementary inclusive analyses have been pursued. In one, candidates for 
ip(2S) — ► n h c — > 7r°(7?7 c ) are selected by choosing events containing an El photon can- 
didate in a range of energies expected for h c — > 7ei^ c , and displaying a peak in M(h c ). This 
method has the advantage that backgrounds to the signal photon and tt° are uncorrelated 
with one another, but it presupposes foreknowledge of the interesting range of M(h c ) values, 
and does not compensate for the broadening of the photon energy spectrum due to h c recoil. 
In a second inclusive method, events are chosen within a given range of M(t] c ) as calculated 
from the energies and relative angle of the tt° and 7ei, and displays a peak in M(h c ). This 
method compensates for the recoil broadening of the 7ei energy spectrum and does not 
presuppose a value of M(h c ). However, since both photon and n° energies are needed to 
calculate M(r] c ), backgrounds are correlated, and some subtraction methods appropriate for 
the first method are not valid for the second. 

Exclusive reconstruction of decay modes of the r] c offers the potential of significant 
background reduction. The following i] c decay modes were studied: KgK^ir 11 , KlK ± 7i T , 
K + K~7t + tt~ , 7T + 7r~7r + 7r~, K + K~ir°, 7r + 7r~r](77), 7t + 7t~7](7t + 7t~7t ). They are summarized in 
Table E together with their branching fractions in r] c decay |7j • In order to reduce the effect 
of the poorly known t] c branching ratios, the ratio of rates of ip(2S) decay to 7r°7?7 c and 7?7 C 
is measured. The normalizing mode has been recently measured at CLEO (3jJ. Its study 
also permits us to construct and verify event selection criteria in rj c reconstruction. 

IV. BACKGROUND SOURCES AND SUPPRESSIONS 

We first describe major backgrounds to the signal, and how they are suppressed, in a 
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TABLE I: Decay modes of n c used in the exclusive analysis and their branching fractions B 



Mode 


B (%) 




1.9 


± 


0.5 


KlK^ 


1.9 


± 


0.5 


K+K--K+-K- 


1.5 


± 


0.6 


7r + 7r~7T + 7T~ 


1.2 


± 


0.3 


K + K-ir° 


1.0 


± 


0.3 




1.3 


± 


0.5 


7r + 7r _ ?7(7r + 7r _ 


-7T°) 0.7 


± 


0.3 


Total 


9.5 ± 1.6 



qualitative manner. Details of background suppression are described in the next section. 
Selection criteria are applied in different ways depending on the nature of the analysis. 

• The transition ip(2S) — > ir + TT~J/?f). Approximately 1/3 of all ip(2S) decay to the 
final state ir + 7T~J/ip |32|]. Subsequent decays of J /if) can generate both soft tt°s (a 
background to the signal for ip(2S) — > 7T°h c ) and hard photons in the vicinity of the 
signal energy E(jei) — 500 MeV for the expected El transition h c — > JeiVc- Thus, all 
analyses to be reported here excluded some range of mass X around M(J/if)) recoiling 
against tt + tt~ in the reaction ip(2S) — > 7t + 7t~X. 

• The transition tp(2S) — > 7T 7T J/if). The decay ip{2S) — > J /if) accounts for about 
1/6 of all if)(2S) decays j3^. In addition to the backgrounds mentioned above for 
charged pion pairs, either of the two neutral pions can be mistaken for that in the 
signal for ip(2S) — > 7i°h c . Thus, in inclusive analyses, a range of masses around 
M(J/ifi) in the spectrum recoiling against the dipion pair in if>(2S) — ► 7r°7r°X was 
excluded. 

• The transition if>(2S) — > 7Xcj — * 77 J/if>- The sum of the product branching ratios 
B(if)(2S) — > lXcj)B{Xcj lJ/4 1 ) exceeds 5% |32 |. This background can be reduced 
by excluding events with a range of masses around M(J/tp) in the spectrum recoiling 
against 77 in 77X. 

• Candidates for 500 Me V El photons which are tt° or rj decay products. A sufficiently 
energetic ir° can give rise to a photon which can be mistaken for the signal El photon 
in h c — > JeiVc- If is possible to suppress such photons by rejecting all candidates which 
can form a candidate 7r° if paired with another photon. A similar rejection of 77 decay 
products also can be applied. 

• Mis-pairings of candidates for tt° decay. In general photons from 7r° decays are iden- 
tified by requiring that their energies and directions lead to a reconstructed 7r° mass 
within about 15 MeV/c 2 of the nominal value of 135 MeV/c 2 . If some other pairing 
gives a better-reconstructed n° mass, the original pairing is discarded and the better 
pairing is adopted. 
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TABLE II: Conditions under which ip{2S) data were acquired for this analysis. Here AE cm denotes 
the center-of-mass energy spread, while / Cdt denotes integrated luminosity measured using the 
reaction e + e~ — ► 77. 



Detector Time AE cm J Cdt N{ip(2S)) 

period (MeV) (pVr 1 (10 6 ) 

CLEO-III 2002-3 L5 2/74 L56 
CLEO-c 2003-4 2.3 2.89 1.52 
Total 5.63 3.08 



V. DATA SAMPLE AND EVENT SELECTION 

The data samples obtained with the CLEO III and CLEO-c configurations are shown in 
Table |TT1 where the number of events was calculated by the method described in (3l[ and 
was estimated to have an uncertainty of ±3%. 

Common features of event selection for all analyses are listed in the following. Several 
other analysis-specific criteria will be described in the corresponding subsections. Selection 
requirements for all analyses are summarized in Table ITTTl 

• Charged particle selection criteria were standard ones used for other CLEO analyses. 
The distance of closest approach of a track with respect to the run-averaged collision 
point was required to be less than 5 cm along the beam line and less than 0.5 cm in the 
direction transverse to the beam. Each track was required to be fitted with a reduced 
X 2 (i.e., per degree of freedom) of less than 20, to give between 50% and 120% of the 
expected number of signals on drift chamber wires, and to make an angle of at least 
21.6° = cos" 1 (0.93) with respect to the beam axis. 

• A photon candidate was defined as a shower which does not match a track within 100 
mrad, is not in a "hot" cell of the electromagnetic calorimeter, and has the transverse 
distribution of energy consistent with an electromagnetic shower. 

• The minimum 7r° photon candidate energy was set at 30 MeV in the barrel and 50 
MeV in the endcaps. 

• In kinematic fitting, photon energies and angles for tt° candidates were adjusted to 
give the exact 7r° mass. This increases precision in the determination of the 7r° energy 
and hence the h c mass, which is computed from Eq. (jSJ) using the nominal values of 
M(i/i(2S)) and M(vr°). 

• Photon candidates for the El transition h c — > ^t\ c were subjected to background 
suppression involving vetoing of candidates which could form a tt°. 

• Neutral pion candidates were tested for the possibility that one of their showers could 
form a neutral pion with some other shower, and were rejected if any other pairing 
was more consistent with a n° mass. 
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• Events were flagged if they were candidates for the processes ip{2S) — > ti + ti J/ip or 
i/j(2S) — > 7r°7r° J '/ip and rejected accordingly. 

• When an empirical parametrization of the background shape was needed, the analy- 
ses employed a convenient parametrization of backgrounds to the ir° recoil spectrum 
known as an ARGUS function [33], appropriate for processes such as ip{2S) — > 7i°h c 
in which there is a kinematic endpoint, equal here to M(ifj(2S)) — M(7r°) = 3551.2 
MeV/c 2 . 

• A large generic Monte Carlo sample of ~ 39 million ip(2S) events permitted the opti- 
mization of signal-to-background ratio by adding an appropriately normalized sample 
of signal Monte Carlo events and choosing event selection criteria to maximize the 
likelihood ratio for fits with and without a resonance signal. 

• The distribution of the photon polar angles in both h c — > ^r\ c and ip(2S) — > ^r\ c 
(relevant to the exclusive analysis) was assumed to be ~ 1+cos 2 9. For the former decay 
this assumption is based on the expectation that the h c retains the ip(2S) polarization 
in the (mainly S-wave) process ip(2S) — > it°h c . 



A. Inclusive analyses 

The event selection criteria for the analysis selecting a range of £'(7ei) are summarized 
in Table 11111 Showers were required to have at least 30 MeV energy if detected in the 
barrel region of the calorimeter and at least 50 MeV if detected in the endcaps. Only 
the ten highest-energy showers and tracks in an event were considered, in order to reduce 
combinatorial background. A maximum of ten neutral pions composed of the ten highest- 
energy showers was considered. 

Neutral pions were reconstructed by requiring that the two-photon invariant mass be in 
the range M 77 = 135 ± 15 MeV/c 2 or within three standard deviations of the peak. (Reso- 
lutions in MeV/c 2 depend on properties of each candidate, such as energy and calorimeter 
location.) 

Selection criteria were guided by maximizing the likelihood ratio for fits to Monte 
Carlo-generated background with and without a simulated signal. In order to reduce 
the abundant background due to photons and charged particles from the decay of J/ip, 
the cascades ip(2S) — > J/tpX were suppressed by excluding candidates for ip{2S) — ► 
(7i + 7i^J/tp, 7T°ir J/'if),'-f'-fJ/ifj) using the criteria in the second column of Table ITTT1 Pho- 
ton candidates for 7 E1 in h c — > JeiVc were rejected if they could form a tt° or r\ (defined, 
respectively, by M 77 = 135 ±15 or 550±25 MeV/c 2 ) when combined with any other photon. 
It was demanded that there be only one photon in the event with energy 503 ± 35 MeV. 

In the complementary analysis selecting a range of M(j] c ) (Table 11111 third column), 
events were chosen corresponding to a slight modification of a previously used criterion [3l| 
for selection of hadronic events at the ip{2S) energy. 1 Background suppression techniques 
were similar in most respects to those of the other inclusive analysis except for the following 
details: 



For 1 < iV c h < 3 (A^ch = number of charged tracks), the maximum energy visible in the calorimeter was 
required to be less than the total center-of-mass energy Eqm, vs. 0.85£"cm in Ref. [il]]. For > 4 the 
criteria were the same as in Ref. |3l| . 
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TABLE III: Comparison of event selection criteria for inclusive and exclusive analyses. 



Property 


Inclusive analysis specifying: 


Exclusive 


or quantity 




A/fin*] ranfp 

1V1 \ l fC) i-aJJ-gC 


analysis 


Initial 


^> 2 cliarp'pH 




Hadronic 


event 


trarks and 


1 i of rhavp'pc] 


selection 


selection 


> 3 showers 


tracks (see text) 


(see text) 


#(7Ei) or 


X -'V /til/ 


M(n^ ± 35 

MVX \'lc ) - 1 - J J 


M(rj c ) ± 50 


M(ry c )range 


503 ± 35 MeV 


MeV/c 2 


MeV/c 2 


Photon 


10 most 


All 


All 


showers 


energetic 






Photon 


Barrel plus 


Barrel 


Barrel plus 


acceptance 


endcaps 


onlv 


endcaps 


No. of 7T° in 


Onp and 


Onp and 


At least 


signal region (a) 


only one 


only one 


one 


7T° rejection 


Reject best- 


Reject all it® 


Reject all ir° 


on TEi 


pull 7T° only 


with pull < 2.5 


with pull < 3 


77 rejection 


M(7Ei7) = 


None 


None 


on TEi 


550 ± 25 MeV/c 2 






|AAf(vr+vr-J/^)| 


< 15 


< 8.4 


< 10 


excluded 


MeV/c 2 


MeV/c 2 


MeV/c 2 


|AA/(vr 7r°J/V')| 


< 40 


< 32 


None 


excluded 


MeV/c 2 


MeV/c 2 




11 J/^ 


Af (all chgd) 


|AAf( 77 J/VOI 


None 


rejection 


> 3050 MeV/c 2 


> 40 MeV/c 2 





(a) Defined as giving M(h c ) = 3526 ± 30 MeV/c 2 

• Photons for 7r° or 7 ei candidates were chosen only in the barrel region of the electro- 
magnetic calorimeter, in an attempt to improve energy resolution. 

• Neutral pion candidates were required to have a 77 mass within 2.5a of the peak, and 
were rejected if any other pairing of photons within this same "pull mass" (normalized 
deviation from the correct mass in units of Gaussian width) provided a better fit to 
the 7T° mass. Partner photons for this rejection were allowed to be either in endcaps 
(E > 50 MeV) or barrel (E > 30 MeV). 

• Candidates for the El transition photon which could form a n° were vetoed j3jj as in 
the £'(7Ei)-selection analysis, rejecting any photon forming a pair with mass less than 
2.5cr from M(tt°) when combined with a photon in endcap regions of the calorimeter 
with at least 50 MeV or barrel regions with at least 30 MeV. However, Monte Carlo 
simulations (to be discussed in Sec. VI) indicated no need to veto rj mesons. 
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B. Exclusive analysis 



The exclusive analysis measures the ratio of the cascade decays i/)(2S) — > n°h c — > -n Q {jrj c ) 
to the direct radiative decays i/j(2S) — > 7?y c by identifying the decay channels listed in Table 
I. To design event selection criteria, 20,000 signal Monte Carlo were generated for each mode 
of the cascade and direct radiative decays. The 39 million generic Monte Carlo ^(2S) decays 
without h c were utilized to study the background to the cascade decay. All reconstructed 
events were required to have no extra tracks and total extra unmatched shower energy less 
than 200 MeV. The basic particle selection criteria, in addition to those mentioned at the 
start of this Section, include the following specific to this analysis: 

• 7T°: Mass less than 3a from nominal value. 

• Kg\ Decay displaced by more than 3a with respect to the run-averaged collision point, 
mass within 10 MeV/c 2 of nominal value 

• v(ll) : Mass within 3a of the nominal rj value 

• ?7(7r + 7r~7r ) : M^+^-^o within 20 MeV/c 2 of the nominal rj mass 

Information from the RICH and dE/dx detectors was combined to distinguish kaons from 
pions when RICH information was available. RICH information was utilized when a track 
was in the RICH fiducial volume with | cos6>| < 0.8, a kaon candidate had momentum at 
least 600 MeV/c, and three or more photons were detected near the predicted ring location. 
A combined "Log-Likelihood" was defined as AL = L(7r) RICH — L(K) K icn + {CdE/dx) 2 ~ 
( a dE/dxY ! > w h ere -^(tOrich is ~2 times the natural logarithm of the RICH likelihood for the 
pion hypothesis, and L(K) rich is for the kaon hypothesis, while a dE / dx is the deviation of 
dE/dx from what is expected for the pion hypothesis normalized to the measurement error 
and a dE / dx is the same for the kaon hypothesis. If RICH information was not available, a 
track was identified as a kaon if \a dE , dx \ < 3 and \CdE/ dx \ < WdE/dx\)- When RICH informa- 
tion was not available and track momentum was above 600 MeV/c, a track was identified 
as a pion if \a dE / dx \ < 3. When RICH information was available or track momentum was 
below 600 MeV/c, charged kaons and pions were well-separated. In the K + K~it + it~ and 
K + K~ir° modes, at least one kaon candidate was required to be identified when K and n 
were well-separated. 

Because the ip(2S) resonance width is only 0.3 MeV, considerably less than the beam 
energy spread, the beam energy was always assumed to be half of M(ip(2S)) when running 
at the vp(2S). 2 In r/ c — > K\K ± Tt^ : , the missing mass should equal the K L nominal mass since 
the K\ is undetected. In this case, a 1C kinematic fit was performed assuming that the 
missing particle has the mass of K\. In all other modes, ip{2S) final decay particles were fully 
reconstructed, and the net 4-momentum of reconstructed charged or neutral tracks should 
equal the 4-momentum of the ip(2S) which is known, permitting 4C kinematic fits. The x 2 
values from the fits indicate how well each reconstructed event matches the kinematics of 
the decay hypothesis. A rather loose requirement of x 2 /d.o.f.< 10 in all modes was imposed. 



2 The crossing angle is around 4 mrad, corresponding to a transverse momentum of about 3686 sin(0.004) = 
15 MeV/c. 
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The r] c signal was fully reconstructed in all the modes except K^K ir^. In K\K ^ y the 
r\ c mass was inferred from the energies of the recoiling El photon and it . 

Generic Monte Carlo studies indicate that photons from 7r°s in ip(2S) — > 7T 7[°J/ip and 
if)(2S) — > JXcj (Xcj —> 7T°X) decays are a large background source to 7ei- A photon 
candidate was vetoed if the absolute value of its best n° pull mass, when combined with all 
other photons of energies greater than 30 MeV, was less than 3. This cut greatly reduced 
the background but also resulted in a 15% efficiency loss according to signal Monte Carlo. 
The net effect on the expected sensitivity to h c was positive. 



Monte Carlo simulations of background and signal were employed in order to optimize 
event selection criteria and to estimate backgrounds to data. The generic Monte Carlo sam- 
ple mentioned earlier was used. Simulations employed hadronization routines embodied in 
JETSET |35j], with its parameters optimized for ip(2S) decays [3l|. The detector simulation 
was based on Geant |36j]. Hadronization of h c decays was emulated using Model 14 of the 
LUND/ JETSET fragmentation algorithm. 

A. Inclusive analyses 

1. Choice of background shapes. 

The .E(7Ei)-range analysis uses the ir° recoil spectrum from the data itself as background, 
without demanding a candidate with E 1 = 503 ± 35 MeV for the El photon. This is feasible 
since the h c contribution is invisible, being at the level of ~ 4 x 10~ 4 . The M(r/ c )-range 
analysis uses generic Monte Carlo background instead, since the selection of an i] c mass 
range in analyzing the data affects the background shape. 

2. Optimization of signal significance. 

Monte Carlo samples were employed to choose ranges of selection providing the highest 
sensitivity to the h c signal, as judged by maximum likelihood for the resonance hypothesis. 
These samples also permitted studies of input/output agreement and statistical variation. 
The optimum event selection criteria determined in these Monte Carlo studies were applied 
to the data. 

In the £'(7Ei)-range analysis, 30,000 signal events were generated for ip(2S) — > 7T°h c — > 
7r°(77/ c ). Assuming B^B h = B(^){2S) -> n°h c ) x B(h c -> jr) c ) = 4.0 x 10~ 4 , 15,600 signal 
events were added to the 39 million generic Monte Carlo sample. The input masses and 
widths were taken as M(h c ) = 3526 MeV/c 2 , T(h c ) = (0.5, 0.9, 1.5) MeV, and M(r] c ) = 2982 
MeV/c 2 , T(i] c ) = 24.8 MeV [30]. In the M(?7 c )-range analysis, 185x 10 3 events were generated 
for ip{2S) — > ir°h c , with a 37.7% branching ratio |23J for the subsequent decay h c — > jrjc. 
The remaining h c decays were taken to have a 56.8% branching ratio to ggg and a 5.5% 
branching ratio to jgg. The mass of h c was assumed to be 3525.3 MeV /c 2 , and the h c width 
was taken to be 1 MeV. The mass of 7] c was chosen as 2981.8 MeV/c 2 |30j. 

The results of the Monte Carlo studies for the _E(7Ei)-range analysis are summarized 
in the second and third columns of Table IIVI Significance levels are obtained as a = 



J — 2 \n(L / L max ) , where L max is the maximum likelihood for the resonance fit, and L 
is the likelihood for the fit with no h c resonance. Selection ranges (summarized in the 
second column of Table IIII|) were chosen to maximize the significance for the Monte Carlo 
sample calculated in this manner. For each effect examined, asterisked values for all other 
parameters were assumed. 



VI. MONTE CARLO SAMPLES 
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TABLE IV: Results of Monte Carlo optimizations using a combined sample of 39 million generic 
i>(2S) events and 15,600 signal events for £?(7Ei)-range analysis. Asterisks show final selection. 





MC 


DATA 




Signif. (a) 


s 2 /B 


Mass, MeV/c 2 


Yield 


B*B h x 10 4 


X 2 /DOF 


Signif. (a) 


Effect of background shapes 


* DATA 
MC 






3524.4±0.7 
3524.6±0.7 


139±41 
146±40 


3.4±1.0 
3.5±1.0 


1.36 
1.59 


3.6 
3.8 


All o 


* the following optimizations were done using background from 


DATA 


Effect of changing range of hard 7 energy, 503±, MeV/c 2 


±30 
* ±35 
±40 
±45 
±50 


16.4 
17.3 
16.1 
16.3 
15.8 


1.01 
1.00 
0.96 
0.90 
0.86 


3524.0±0.7 
3524.4±0.7 
3524.4±0.6 
3524.8±0.8 
3524.8±0.9 


120±38 
139±41 
145±43 
134±45 
132±47 


3.1±0.9 
3.4±1.0 
3.4±1.0 
3.0±1.0 
2.9±1.0 


1.19 
1.36 
1.28 
1.24 
1.26 


3.3 
3.6 
3.5 
3.1 
2.9 


Effect of c 


ranging mass range for tt + it J/tp rejection, MeV/c 2 


±6 
±10 
* ±15 
±20 


17.2 
17.3 
17.3 
17.1 


0.97 
0.99 
1.00 
1.00 


3524.4 ± 0.6 

3524.3 ± 0.6 

3524.4 ± 0.7 
3524.2 ± 0.7 


158 ± 43 
156 ± 42 
139 ± 41 
132 ± 40 


3.7 ±1.0 
3.7 ±1.0 
3.4 ± 1.0 
3.3 ±1.0 


1.28 
1.36 
1.36 
1.38 


3.9 
3.9 
3.6 
3.4 


Effect of changing mass range : 


"or 7r°7r°J/V ; rejection, MeV/c 2 


±20 
±30 
* ±40 
±50 


17.2 
17.2 
17.3 
17.3 


0.99 
1.00 
1.00 
1.00 


3524.3 ± 0.8 
3524.5 ± 0.8 

3524.4 ± 0.7 
3524.4 ± 0.7 


140 ± 42 
134 ± 41 
139 ± 41 
147 ± 41 


3.3 ±1.0 
3.2 ±1.0 

3.4 ±1.0 
3.6 ± 1.0 


1.45 
1.30 
1.36 
1.30 


3.4 
3.4 
3.6 
3.8 




Effect of number of 7r°s in the signal region 


* = 1 
> 1 


17.3 
17.2 


1.00 
0.95 


3524.4 ± 0.7 
3524.8 ± 0.9 


139 ± 41 
122 ± 42 


3.4 ±1.0 
2.9 ± 1.0 


1.36 
1.04 


3.6 
3.0 


Effect of endcap 7s in signal 7r°s 


* with 
without 


17.3 
16.0 


1.00 
0.91 


3524.4 ± 0.7 
3524.8 ± 0.7 


139 ± 41 
123 ± 37 


3.4 ± 1.0 
3.4 ± 1.0 


1.36 
1.16 


3.6 
3.5 


Effect of 77 suppression on El photon 


* with 
without 


17.3 
15.8 


1.00 
0.91 


3524.4 ± 0.7 
3524.6 ± 0.8 


139 ± 41 
135 ± 45 


3.4 ± 1.0 
3.0 ±1.0 


1.36 
1.21 


3.6 
3.1 


El 


feet oi 


" i/)(2S) — > 7x1,2 — ► 11 J W suppression 


* without 
with 


17.3 
17.0 


1.00 
1.02 


3524.4 ± 0.7 
3524.6 ± 0.7 


141 ± 41 
137 ± 40 


3.4 ± 1.0 
3.4 ± 1.0 


1.36 
1.21 


3.6 
3.6 
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Table IV, continued 






DATA 


'no MC entries) 






Mass, MeV/c 2 


yield 


$A x io 4 


X7DOF 


signif. (a) 




Effect of chan 


ging total width of h c , MeV 




0.5 


3524.3 ±0.7 


132 ± 38 


3.2 ±0.9 


1.36 


3.6 


* 0.9 


3524.4 ±0.7 


139 ±41 


3.4 ± 1.0 


1.36 


3.6 


1.5 


3524.5 ±0.7 


149 ± 44 


3.6 ± 1.1 


1.39 


3.6 




Effect of changing tt° 


resolution widths 




* MC 


3524.4 ±0.7 


139 ±41 


3.4 ± 1.0 


1.36 


3.6 


MC-25% 


3524.3 ±0.6 


131 ±38 


3.2 ±0.9 


1.35 


3.6 


MC+25% 


3524.5 ±0.7 


149 ± 45 


3.6 ± 1.1 


1.39 


3.6 


Effect of binning 


* 2 MeV/c 2 


3524.4 ±0.7 


139 ±41 


3.4 ± 1.0 


1.36 


3.6 


1 MeV/c 2 


3524.5 ±0.6 


137 ±41 


3.3 ± 1.0 


1.16 


3.5 


Effect of changing fii 


range, MeV/c 2 


* 3496-3552 


3524.4 ±0.7 


139 ±41 


3.4 ± 1.0 


1.36 


3.6 


3500-3540 


3524.4 ±0.7 


139 ± 42 


3.4 ± 1.0 


0.96 


3.5 


CLEO III VERSUS CLEO-c 


CLEO III 


3523.8 ±0.7 


94 ±30 


4.5 ± 1.4 


0.96 


3.3 


CLEO-c 


3526.1 ± 1.5 


56 ±28 


2.8 ± 1.4 


1.55 


2.1 



These choices were found to lead to the same output from the if)(2S) generic Monte Carlo 
sample as the input: 





Input 


Output 


M(h c ) (MeV/c 2 ) 
B^B h x 10 4 


3526.0 
4.0 


3525.9 ±0.1 
4.1 ±0.3 



The above choices were based on maximum likelihood in 22 variations with no contact 
with the experimental data, i.e., by "blind" analysis. The best choices indeed are mirrored in 
the data. Table HVI therefore lists for the data the values of the likelihood-based significance 
for all 22 variations examined in the Monte Carlo sample. It is interesting to note that these 
choices do lead to higher significance values in most cases, although, as is to be expected, 
because of the factor ~ 13 smaller statistics in the data, both the significance level and their 
variations are smaller than those in the Monte Carlo sample by a factor close to \/l3. 

The 7T° recoil mass distribution for the Monte Carlo sample in the £'(7Ei)-range analysis 
is shown in Figure 121 It was fitted using the sum of two Gaussians with widths fixed to 
values determined by the signal Monte Carlo sample. The background was fitted using a 
histogram of the ir° recoil distribution from the generic Monte Carlo as described above. 
The dashed line shows the contribution of background without signal. 

In the M(^ c )-range analysis, widths in M(rj c ) were determined by fits using a Gaussian 
plus a low-order polynomial, while fits to M(h c ) used a Breit-Wigner resonance function 
with T = 1 MeV convolved with two Gaussians, a quadratic polynomial constrained to 
vanish at the kinematic endpoint, and an ARGUS background function. 

The best range of 7] c masses for optimizing signal significance was determined via Monte 
Carlo studies using a likelihood ratio criterion. Five M(r) c ) windows 2940-3020, 2945-3015, 
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FIG. 2: Spectrum of masses (in GeV/c 2 ) recoiling against 7r° in a sample of 39 million generic 
Monte Carlo events plus 15600 signal Monte Carlo events (i?(7Ei)-range inclusive analysis). The 
solid histogram illustrates the fit described in the text. 



2950-3010, 2955-3005, and 2960-3000 MeV/c 2 were considered. Upper and lower bounds 
were chosen symmetrically with respect to M(r] c ) ~ 2980 MeV/c 2 . Detection of the correct 
candidate for the El photon but assignment of a background n° with the wrong energy as a 
signal 7T° candidate can introduce a potential bias on M{h c ) in the presence of asymmetric 
M(r] c ) limits. 

Selecting events within the above M(r] c ) windows, fits were performed for 3496 MeV/c 2 < 
M(h c ) < 3551.2 MeV/c 2 to the generated h c mass distributions. The signal Monte Carlo 
was generated using a flat angular distribution for the El photon. A correction to the 
efficiency was performed for the expected form W^(cos6>) ~ 1 + cos 2 9 with respect to the 
beam axis. The ratio of the two efficiencies when integrating to a maximum |cos# max | is 
ReS = (1/4) (3 + cos 2 $ max ). For |cos# max | = 0.804, corresponding to the outermost ring 
of the barrel calorimeter used in this analysis, the correction factor is i? e fr = 0.912. The 
efficiencies were corrected for R c q. 

After fits to the signal Monte Carlo yielded the parameters of its Breit-Wigner plus 
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FIG. 3: Generic Monte Carlo M{h c ) distribution (M(r/ c )-range inclusive analysis) for simulated 
ip(2S) data of 39 million events with a signal of 69.7 x 10 3 h c decays corresponding to 15.8 x 10 3 
events of h c — ► 7r/ c for 2945 MeV/c 2 < M(rj c ) < 3015 MeV/c 2 . The generated masses were 
[M(r/ C ), M(h c )] = (2981.8, 3525.3) MeV/c 2 . The signal was emulated using a pair of Gaussians and 
a Breit-Wigner with T = 1 MeV. The dashed line shows the contribution of background. 



Gaussian functions, the generic Monte Carlo distribution was combined with a weighted 
signal distribution to emulate a combined branching ratio for the decay ip(2S) — > it h c 
followed by h c — > JeiVc of B^Bh = 4 x 10~ 4 . The resulting distribution was fitted both 
with (generic + weighted signal), and with generic background alone, yielding a ratio of 
likelihoods. 

This process resulted in an optimum range of 2945 MeV/ c 2 < M(r] c ) < 3015 MeV/c 2 . The 
corresponding M(h c ) distribution is shown in Fig. El For any wider M(r) c ) range, the photons 
from the transition h c — > 7^ c become contaminated with contributions of Doppler-broadened 
photons from the El transition Xc2(3556) — > ^yJ/ip- Backgrounds from this transition and 



17 



TABLE V: Fits to simulated signal and background using a Breit-Wigner signal function convolved 
with a double Gaussian and a generic Monte Carlo background (M(r/ c )-range analysis). Branching 
ratios include an efficency factor R c s = 0.912 for the 1 + cos 2 8 distribution of the El photon. The 
nominal M{g c ) range is labeled by an asterisk (*). 





M(r? c ) range (MeV/c 2 ) 


2940-3020 


*2945-3015 


2950-3010 


2955-3005 


2960-3000 


M(/i c )(MeV/c 2 ) 
Significance a 
Efficiency (%) 
B^Bh x 10" 4 


3525.24±0.16 
17.08 
15.3 
4.07 ±0.25 


3525.23±0.16 
17.30 
14.6 
4.07 ±0.25 


3525.22±0.17 
17.20 
13.5 
4.07 ±0.25 


3525.21±0.17 
17.05 
12.2 
4.07 ±0.25 


3525.18±0.18 
16.45 
10.6 
4.07 ±0.26 



others rise steeply as the upper limit on M(r] c ) is increased above 3020 MeV/c 2 . 

Fits to simulated signal and background in the M(?7 c )-range analysis are compared in 
Table fVl The r] c mass range 2945-3015 MeV/c 2 gives the greatest signal significance for an 
h c of mass 3525.3 MeV/c 2 produced with B{ip{2S) -> 7i°h c )B{h c -> 777J = 4 x 10" 4 . The 
extracted values of M(h c ) are about 0.1 MeV/c 2 below the input. This feature is included 
in the estimate of systematic errors. The maximum significance of 17.3a scales to 4.8a for 
a sample of 3.08 x 10 6 events. 

3. Variations in output parameters. 

In the generic Monte Carlo sample, for all the 22 variations of the £'(7Ei)-range analysis 
listed in Table ITVl the change in output M(h c ) and B^Bh were found to be AM(h c ) < 0.1 
MeV/c 2 , and A(B^Bh) < 0.2 x 10~ 4 , i.e., within the statistical errors assigned by the output. 
To see the level of statistical variations in Monte Carlo samples as small as the data (i.e., ~ 3 
million '0(25')), the total sample of 39 million ip(2S) decays was split into 13 independent 
samples, each of 3 million ip(2S). Table IVT1 summarizes results of the analysis for the choices 
of the final selection and for variations of these choices. For the final selection the limits of 
variation were found to be AM = (-0.4, ±0.3) MeV/c 2 and A(B^B h ) = (-1.1, ±1.4) x 10" 4 . 
For B^Bh the effect of variations from the final selection is within the range observed for 
the final selection. There may be some evidence of larger than expected variation when one 
changes AJ5(te;l) to ±50 MeV, and when one includes more than one signal tt° candidate. 
A choice of AE(jei) — 50 MeV begins to accept photons on the high-energy tail of the 
transition Xc2 ~^ lJ when detector resolution and recoil effects are taken into account. 

Because the Monte Carlo signal sample was generated with an assumed M(h c ) = 3526 
MeV/c 2 , or Efai) = 503 MeV, it is prudent to examine what bias is introduced in M(h c ) 
and B^Bh if the true M(h c ) were to differ from 3526 MeV/c 2 . The resulting variation in 
efficiency was found to be less than 2.5% for M(h c ) = 3526 ± 14 MeV/c 2 . 

The corresponding variations in the M(7/ c )-range analysis were explored by again forming 
13 samples of ~ 3 million generic ip(2S) Monte Carlo and adding 13 samples of 3135 signal 
Monte Carlo events with B(h c — > -yr] c ) = 37.7%, B(h c — > ggg) = 56.8%, and B(h c — > -ygg) = 
5.5%. This permitted simulation of a combined branching ratio B^Bh = 4 x 10~ 4 . Fits were 
performed using the same functions used in fitting data. The results are shown in Table IVTT1 
Deviations from the mean were found to be of the expected magnitude for data samples of 
this size. 
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TABLE VI: Results for M(h c ) and B^Bh from trial experiments with 13 independent Monte Carlo 
samples of 3 million ip(2S) each [£'(7Ei)-range analysis]. The inputs were M(h c ) = 3526.0 MeV/c 2 
and B^B h = 4.0 x 10" 4 . The full Monte Carlo sample yielded M(h c ) = 3526.1 ± 0.1 MeV/c 2 and 
B^B h = 4.1±0.3 x 10~ 4 . Variations from the final selection resulted in AM(h c ) < 0.1 MeV/c 2 and 
A(B 1 pB h ) < 0.2 x 10~ 4 for this large sample. The second column lists AM(h c ) = M(h c ) - 3526 
MeV/c 2 or A(B i pBf l ) = (B^Bh) — 4.0 x 10~ 4 for the final selection. The following columns list 
AM{h c ) or A(B^pBf l ) for the specified variations from the final selection. The statistical error on 
all output masses was ±0.5 to ±0.6 MeV/c 2 and on all output B^Bh was 1.0 x 10~ 4 . 





AM(h c ) - MeV/c 2 


AM(h c ) - MeV/c 2 with variations from final 




Final selection 


AE 1 ± 50 MeV 


> lvr 


No endcap 


No n supp. 


MC 


-0.4/+0.3 


-1.8/+0.7 


-2.1/+0.2 


-0.3/+0.2 


-0.4/+0.3 


Data 




±0.4 


±0.4 


±0.4 


±0.2 




A{B^B h x 10 4 ) 


A(B i f,Bh x 10 4 ) - with variations from final 




Final selection 


A£ 7 ± 50 MeV 


> lvr 


No endcap 


No n supp. 


MC 


-1.1/+1.4 


-1.2/+0.7 


-0.3/+1.1 


-0.5/+0.3 


-1.1/+0.3 


Data 




-0.5 


-0.5 


±0.0 


-0.4 



TABLE VII: Results for M{h c ) and B^B^ from trial experiments with 13 independent Monte Carlo 
samples of 3 million ip(2S) each [M(ry c )-range analysis]. The inputs were M(h c ) = 3525.3 MeV/c 2 
and B^B h = 4.0 x 10" 4 . The full Monte Carlo sample yielded M(h c ) = 3525.33 ± 0.18 MeV/c 2 
and B^Bh = 3.9 ± 0.3 x 10" 4 . The second column lists AM(h c ) = M(h c ) - 3525.3 MeV/c 2 or 
A(B^Bh) = (B^pBh) — 4.0 x 10~ 4 for the final selection. The following columns list AM(h c ) or 
A(B^Bh) for variations from the final selection. 





AM(h c ) (MeV/c 2 ) 
Final selection 


AM(h c ) (MeV/c 2 ) with variations from final 


AM(i] c ) ± 40 MeV 


AM(t) c ) ± 20 MeV 


>l7T U 


w/endcap 


w/77 supp. 


MC 
Data 


-0.5/+0.3 


-0.5/+0.3 
±0.4 


-0.5/+0.4 
-0.3 


-0.4/+0.3 
±0.0 


-0.4/+0.4 
±0.5 


-0.4/+0.5 
-0.4 




A(B^B h x 10 4 ) 
Final selection 


A(B^Bh x 10 4 ) - with variations from final 


AM(i] c ) ± 40 MeV 


AM(t) c ) ± 20 MeV 


> lvr 


w/endcap 


w/77 supp. 


MC 
Data 


-0.7/+0.5 


-0.7/+0.4 
-0.6 


-0.6/+0.4 
±0.1 


-0.6/+0.3 
-0.3 


-0.5/+0.4 
-1.5 


-0.5/+0.6 
±0.1 



4. Quality of generic Monte Carlo simulation. 

Because the CLEO generic Monte Carlo is used to determine optimum selection criteria 
for energy ranges and binary choices, one must quantify its level of agreement with data in 
emulating the M(h c ) spectrum. The EvtGen |24| generator is combined with a JETSET 
[35l | version tuned to match the relevant low-energy regime j^]. For photon energies below 
450 MeV and pion momenta below 550 MeV/c, the data and Monte Carlo agree within ±5%. 
Above these values the ratio of data to Monte Carlo falls below 95%, rising again from ~ 90% 
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TABLE VIII: Binary choices of selection and criteria (M(r/ C ) analysis). Asterisks denote nominal 
choices. 



Ill 

mTk 


Range 


V 

supp 


MC 
a 


Signal 


Mass (MeV/c 2 ) 


Evts. in pk. 


B (KT 4 ) 


*77 


*M( Vc ) 


*No 


17.3 


3525.3 ± 0.6 


159 ± 41 


3.5 ±0.9 


*77 


*M(i lc ) 


Yes 


16.9 


3524.9 ± 0.6 


132 ± 35 


3.6 ± 1.0 


*77 


£(7El) 


*No 


16.7 


3525.3 ±0.7 


161 ± 44 


3.4 ±0.9 


*77 


#(tei) 


Yes 


16.3 


3524.8 ± 0.6 


134 ± 37 


3.6 ± 1.0 


mTk 


*M(i 1c ) 


*No 


17.3 


3525.1 ±0.6 


152 ± 42 


3.3 ±0.9 


mTk 


*M( Vc ) 


Yes 


16.9 


3524.7 ± 0.6 


134 ± 36 


3.6 ± 1.0 


mTk 


E(jei) 


*No 


16.6 


3525.1 ±0.7 


145 ± 41 


3.1 ±0.9 


mTk 


^(7El) 


Yes 


16.2 


3524.7 ± 0.5 


136 ± 38 


3.6 ± 1.0 



above E 1 = 600 MeV and from ~ 85% above p(tt°) = 950 MeV/c. For low energy photons 
in the slow 7r° from ip(2S) — > 7r°h c , the generic Monte Carlo is satisfactory, but its use over 
extended ranges of energy and momenta, as required in determining background shapes, 
may not be so. This provides a motivation for basing the background shapes on the data, 
i.e., the ir° recoil spectrum without requiring E 1 = 503 ± 35 MeV, instead of the 7r° recoil 
spectrum from the generic Monte Carlo. 

5. Choices in M(j] c ) analysis. 

In the £"(7ei) analysis, electromagnetic cascades involving El transitions to and from 
intermediate Xc states were suppressed by excluding events with the effective mass of charged 
tracks exceeding 3050 MeV/c 2 ("mTk" criterion). In the M(j] c ) analysis, the mass recoiling 
against 77 was reconstructed directly ("77" criterion), and events with a recoil mass within 
±40 MeV/c 2 of M(J/ip) were excluded. 

In the M(rj c ) analysis, which does not use endcap photons and does not restrict photons 
in 7T° candidates to the ten most energetic showers, an advantage in Monte Carlo significance 
by about 0.6a appears when the M(j] c ) range rather than the £?(7ei) range is selected. 

In the .E(7ei) analysis, Monte Carlo likelihood ratios favor suppressing 7E1 candidates 
which can form an rj when paired with other photons. In the M(rj c ) analysis, which uses a 
larger pool of photon candidates for possible pairings, such a suppression entails a loss of 
efficiency for signal detection, leading to decreased significance in Monte Carlo by 0.4cr. The 
M(t) c ) analysis consequently does not adopt this suppression. 

The above three criteria were compared in a binary manner, leading to the results shown 
in Table I Villi The effects of each variation are largely independent of each other when 
measured by change in significance. The first row was chosen over the fifth in the M(j] c ) 
analysis on the basis of a very slight excess in Monte Carlo (MC) significance a; differences 
in resulting mass and branching ratio are within statistics. 

6. Dependence on branching ratio B(h c — > 7ei^ c ) and M(h c ) in signal Monte Carlo. 

In the E(jei) analysis, Monte Carlo simulations were performed by assuming Bh = 
B(h c — > IeiVc) = 100% rather than the value of 37.7% [23j used in the M(rj c ) analysis. 
Moreover, slightly different values of M(h c ) for the signal Monte Carlo were used in the two 
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analyses. The results of changing just Bh or both Bh and M(h c ) in the signal Monte Carlo 
were studied for the M(rj c ) analysis. Several features were notable in this comparison. 

(1) The maximum signal likelihoods in Monte Carlo were less for the choice of Bh = 
100%: (15.5,16. l)a for M(h c ) = (3525.3,3526.0) MeV/c 2 versus 17.3a for B h = 37.7% and 
M(h c ) = 3525.3 MeV/c 2 . (2) For the same M(i] c ) range, the values of M(h c ) in data were 
stable under variation of Bh or input M(h c ), while the extracted values of B^Bh rose by about 
0.4 x 10~ 4 when B h = 100% was taken in the signal Monte Carlo. (3) When B h = 100%, 
the maximum signal likelihood in Monte Carlo still favored no 7] suppression applied to the 
El photon, but to a lesser extent. 

Because the variations in M(h c ) and B^Bh observed under the above changes were ascrib- 
able to the signal fitting hypothesis rather than to the data themselves, they were included 
in estimates of systematic error, giving AM(h c ) = —0.1 MeV/c 2 and AB^Bh = +0.4 x 10~ 4 . 

7. Asymmetric M(t] c ) selection windows. 

The 7] c mass windows were chosen symmetric about 2980 MeV/c 2 in the M(rj c ) analysis 
to avoid M(h c ) spectrum distortions if an El photon of the correct energy were paired with a 
random pion not associated with the transition ip(2S) — > 7r°h c . Slightly higher Monte Carlo 
significance (17.5cr versus nominal 17.3<r) occurs with the asymmetric window 2955-3015 
MeV/c 2 (versus nominal 2945-3015 MeV/c 2 ). On the other hand, the signal significance 
in data peaks for the asymmetric window 2945-3005 MeV/c 2 at 4.6a (versus 4.0<r for the 
nominal window), and the value of M(h c ) obtained from the data is 0.4 MeV/c 2 lower. 
This behavior is consistent with the lower r] c masses observed in a recent analysis of ip{2S) 
radiative decays [3l[ and in the exclusive analysis reported below. 

B. Exclusive analysis 

The signal Monte Carlo indicates that the reconstructed (or recoil) i] c mass and width are 
mode dependent because of the different final decay particles. The value of M(r] c ) calculated 
after kinematic fitting was required to be within 50 MeV/c 2 of the nominal mass. Monte 
Carlo events indicate that this is more than 80% efficient. The width of the reconstructed 
i] c mass distribution depends on both the detector resolution and the intrinsic width, T(rj c ). 
The latter has not been well measured 0, IHj, and the former is decay-mode dependent. 
Because the requirement that M(rj c ) be within 50 MeV/c 2 of its nominal value is loose, 
the systematic uncertainty of the efficiency due to this requirement is minimal, however. 
Measuring the ratio of branching ratios for cascade decay and direct radiative decay reduces 
this systematic uncertainty further. In addition to the other criteria in Table III, this analysis 
takes the 7r° pull mass limit for signal selection and 7i~° suppression to be 3, and the reduced 
X 2 for kinematic constraints to be less than 10. The direct radiative decay ip{2S) — > jr) c is 
studied in the same rj c decay modes, using similar event selection criteria except that the 
M(r] c ) and signal tt° selection criteria are dropped, and the r\ c yield is determined from the 
fit to the 7 recoil mass spectrum. 
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FIG. 4: M(h c ) distribution from recoil 7r° for the CLEO III + CLEO-c data set corresponding to 
the final event selection in inclusive analysis based on selecting a range of £'(7ei). The dashed 
line denotes the background function. The x 2 P er degree of freedom for the fit including peak and 
background is 34.1/25 = 1.36, as noted in Table HVI The corresponding confidence level is 10.5%. 



VII. THE SIGNAL IN THE DATA 



A. Inclusive analyses 



Figure |U shows the spectrum of recoils against 7r° for the data in Table ITT1 with the event 
selection criteria determined to optimize the signal sensitivity in the £'(7ei) analysis. These 
data were fitted with background as determined in Sec. VI plus a Breit-Wigner resonance 
of width 0.9 MeV. The background used was the 7r° recoil spectrum without the cut on 
E('~fEi)- The peak shape consisted of the Breit-Wigner width convolved with an instrumental 
resolution function, determined from the signal Monte Carlo simulation, which itself was 
fitted with a double Gaussian. The efficiency for the final event selection was determined to 
be e = 13.4%. The results are: 



N(evts) = 139 ±41, significance = 3.6a 
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FIG. 5: M(h c ) distribution from recoil tt° for 2945 MeV/c 2 < M(r] c ) < 3015 MeV/c 2 , fitted over 
the range 3496 MeV/c 2 < M(h c ) < 3551.2 MeV/c 2 [analysis selecting range of M(rf c )]. The curve 
denotes the background function based on generic Monte Carlo plus a signal as described in Sec. 
VI B. The dashed line shows the contribution of background alone. The peak contains 159 ± 41 
events. The confidence level of the fit to signal + background was 34%, corresponding to x 2 = 55.6 
for 52 degrees of freedom. 



M(h c ) = 3524.4 ± 0.7 MeV/c 2 
B^B h = B{^{2S) - n°h c ) x B{h c 



7Vc 



(3.4±1.0)xl0~ 4 . 



When selecting a range of M(?7 c )from Monte Carlo, choosing events in the interval 2945- 
3015 MeV/c 2 gave the greatest signal significance, and hence this interval was used for 
further analysis. For the data the significance is slightly greater for a narrower range of 
M(rj c ), as shown in Table HXl The resulting h c mass spectrum is shown in Fig. |S| The 
results are: 



N(evts) = 159 ±41, significance = 4.0<r 
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TABLE IX: Same as Table |V| for fits to CLEO-III and CLEO-c tp{2S) data [M(r] c ) analysis]. 





M(r] c ) range (MeV/c 2 ) 


2940-3020 


*2945-3015 


2950-3010 


2955-3005 


2960-3000 


M(/i c )(MeV/c 2 ) 
Signif. a 
BfBh x 10 4 


3525.67±0.85 
3.24 
2.86 ±0.91 


3525.26±0.60 
4.03 
3.53 ±0.91 


3525.08±0.55 
4.27 
3.76 ± 0.92 


3525.06±0.57 
4.22 
3.76 ±0.93 


3524.97±0.58 
3.97 
3.65 ±0.97 



TABLE X: M(h c ) and combined branching ratio B^B\ t for separate CLEO-III and CLEO-c data 
samples \M{j] c ) analysis, range 2945-3015 MeV/c 2 ]. 



Data 
sample 


Mass 
(MeV/c 2 ) 


Events 
in peak 


Branching 
ratio (10 -4 ) 


CLEO-III 
CLEO-c 


3524.1 ± 1.0 
3526.6 ±0.8 


86 ±29 
93 ±29 


3.8 ± 1.3 
4.2 ± 1.3 



• M(h c ) = 3525.3 ± 0.6 MeV/c 2 

• B^B h = (3.5 ± 0.9) x 10~ 4 . 

The CLEO-III and CLEO-c data were fitted separately. Results are shown in Table ITVl for the 
E^ei) analysis and Table IXI for the M(r] c ) analysis. The relative weights of the two samples 
[with values of M(h c ) differing by about 2 MeV/c 2 ] differ between the two analyses, with the 
-S(7ei) analysis finding fewer signal events in the CLEO-c sample while the M(r] c ) analysis 
finds approximately equal signals in the CLEO III and CLEO-c samples. This accounts 
for the major part of the difference between M(h c ) values in the combined samples. No 
such difference was found in Monte Carlo simulations of CLEO-c data, indicating that the 
observed difference is purely statistical. 

The angular distribution of the 7ei photon in the inclusive analysis was obtained by 
fitting separately the h c peak in the angular ranges 0.0 < | cos#| < 0.3, 0.3 < | cos^l < 0.6, 
and 0.6 < | cos^l < 0.9. The results are presented in Fig. |3 A 1 + cos 2 6 distribution, as 
expected for an El transition from a spin 1 state, gives a satisfactory fit, with \ 2 — 1-7 f° r 
2 degrees of freedom. The angular distribution for the background, obtained in the same 
way as for the fit to the signal, corresponds to the dotted histogram in Fig. HJ and is flat as 
expected. 

B. Exclusive analysis 

There are several ways to search for an h c signal in exclusive modes. One may observe 
enhancements in the photon energy spectrum from h c — > 7%, the reconstructed h c mass 
spectrum, or the recoil 7r° energy spectrum. The photon energy resolution a(E)/E is 2.1% 
to 3.8% for a photon of energy around 500 MeV, depending on whether it is in the barrel or 
endcap Csl calorimeter. The signal photon energy also has a spread because of the intrinsic 
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FIG. 6: Angular distribution of the photons with E~ = 503 ± 35 MeV from the inclusive analysis. 
Solid points denote yield of photons from h c — ► ^yrj c , while open circles denote background photons. 
The curve shows the fit of the h c — > 7% points with a 1 + cos 2 6 distribution. The background 
photons are seen to be isotropically distributed. Scales for the three plots are arbitrary. 



width of 7] c . The reconstructed h c mass calculated from the 4-momenta of the rj c and the 
transition photon also has poor resolution, and depends on r\ c decay modes. In the signal 
Monte Carlo, both the photon energy resolution and reconstructed h c mass resolution are 
larger than 15 MeV in all modes used. The recoil tt° (from ip(2S) — > 7T°h c ) has much better 
energy resolution because of the 7r° mass constraint fit employed in the 7r° reconstruction 
algorithm, as mentioned previously. The M(h c ) spectrum recoiling against a tt° is also 
independent of r] c decay modes, so one can fit the h c signal with the same signal shape when 
signals from different modes are added together. 

After all the selection criteria except for M(t] c ) are imposed, there is a clear cluster of 
events in the plot of r] c candidate mass versus 7r° recoil mass, shown in Fig. [7| Properties 
of the nineteen events in the M(rj c ) band between the dotted lines and with M(h c ) between 
3516 and 3530 MeV/c 2 are summarized in Table IXT1 

There is a highly populated band at the J j ip mass in Fig. [7| Monte Carlo studies indicate 
that most of these events are from 7r°7r°J/'?/> and 7Xcj(</ = 0, 1,2). When one soft photon 
from a 7r° of n°ir J/ip is missing, neither the beam energy constraint nor tt° suppression can 
remove this background, but t] c mass selection is powerful in rejecting such events. Once 
this selection is imposed, corresponding to the range M(j] c ) = 2982 ± 50 MeV/c 2 in Fig. |Hl 
a clearer h c signal appears in the ir° recoil mass spectrum around 3525 MeV/c 2 (Fig. [SJ). 
The distribution was fitted using an unbinned maximum likelihood method and ARGUS 
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FIG. 7: Scatter plot of the reconstructed rj c mass versus the h c candidate mass obtained from tt° 
recoil in data for the exclusive analysis. The horizontal band near M(J/0) = 3097 MeV/c 2 and the 



IXcO, 



diagonal band at larger rj c candidate mass correspond to V>(2S) — > 7r tt°J/iJ) and 4>(2S) 
respectively. The dashed lines denote the region M(r) c ) = 2982 ±50 MeV/c 2 . In this band a cluster 
of events is visible around M(h c ) = 3524 MeV/c 2 . 



background function to obtain the yield and the mass of the observed h c signal. The double 
Gaussian signal shape is obtained from signal Monte Carlo in which the dominant narrower 
Gaussian width is 3.2 MeV/c 2 . The unbinned maximum likelihood fit yields 17.5 ± 4.5 h c 
candidates with mass at 3523.6 ±0.9 MeV/c 2 . The significance of the signal calculated from 
the difference in the likelihood with and without the signal contribution is 6.1a. 

A clear rj c signal also is observed in mass recoiling against the photon in the study of the 
radiative decay ip(2S) — > jrjc. This confirms the appropriateness and effectiveness of the 
event selection criteria. The recoil mass resolution is identical for all modes, and independent 
of track momentum resolution. The signal shape function, a Breit-Wigner function convolved 
with a double Gaussian, is obtained from signal Monte Carlo. The width of the Breit-Wigner 
function represents the rj c intrinsic width. The detector resolution, represented by a double 
Gaussian, was obtained by fitting the distribution of the difference between the generated 
and reconstructed 7] c candidate masses. 

A total of 220 ± 22 events in all seven modes was observed (Fig. [TUJ). The ratio of 
the branching ratios B for the cascade (0(25) — > n°h c — > 7r°(7?7 c )) and direct radiative 
('0(25') — > 7^ c ) decays in each mode is shown in Table IXTT1 To calculate the resulting event- 
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TABLE XI: List of exclusive event candidates. 



Mode 


M(h c ) 


El 


M(rj c ) (MeV/c 2 ) 




(MeV/c 2 ) 


(MeV) Reconstructed Recoil 




3524.3 


475.0 


3018.7 


3012.0 




3529.3 


496.4 


2995.3 


2991.9 


KlK^ 


3521.7 


513.4 




2964.2 




3521.5 


541.2 




2930.8 




3517.7 


463.2 




3019.2 




3523.5 


486.1 




2998.3 




3525.0 


499.9 


2989.2 


2983.4 




3524.3 


474.5 


2978.8 


3012.7 




3526.7 


507.1 


2989.5 


2976.8 


7T + 7r~7r + 7r~ 


3527.2 


494.1 


2983.3 


2992.6 




3520.4 


475.9 


2975.3 


3007.1 




3523.0 


471.6 


2987.5 


3014.8 




3530.9 


523.0 


2956.5 


2962.0 




3519.2 


498.7 


2992.6 


2979.0 




3519.8 


463.2 


3009.1 


3021.3 




3524.0 


473.8 


3007.6 


3013.2 




3524.8 


517.5 


2972.5 


2962.4 


K+K--K 


3525.4 


497.7 


2976.1 


2986.5 


tt + n~ rj(ir + tt~ n®) 


3521.1 


414.4 


3013.0 


3078.8 



weighted average ratio B^Bh/B(dir) = B(ip(2S) — > n h c — > ir ('yrj c ))/B(ip(2S) — > jr] c ), one 
may write the observed number N(X, h c ) of rj c decays via ip(2S) — > 7r°h c — > 7r ( , yr) c ) and the 
observed number N(X, dir) via ip{2S) — > •yr] c to an rj c channel X with B[r\ c — > X) = B(X) 
respectively as 

N(X, h c ) = B i ,B h B(X)N(^(2S))e(X, h c ) , N(X, dir) = B (dir) B (X)N(^j (2 S))e(X, dir), 

(7) 

where e(X, dir) and e(X, h c ) are efficiencies for mode X for direct and cascade decays (Table 
IXHj) . One then finds 

J^TT - - 0-^8 ± 0.049 (stat), (8) 

B(dir) Y,xN{X,drx)'Y,xe{X,dYx)B{X) K h K> 

where Ex N(X, h c ) = 17.5 ± 4.5 and £ x dir) = 220 ± 22. 
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TABLE XII: Efficiencies and yields of direct radiative decay (tp(2S) — > jrj c ) and cascade decay 
(ip(2S) — > vr°/i c — > vr°(777 c )) in exclusive analysis, and ratio of branching ratios, for each mode. 



Mode 


direct radiative decay 


cascade decay 


S(cascade) / 
i3(direct) 


Eff(%) 


Yield 


Eff(%) 


Yield 




12.7 


35.5 ±7.6 


5.6 


1.9±1.4 


0.116±0.090 




32.6 


74.0±12.0 


15.3 


3.1±2.1 


0.081±0.057 


K+K-TT+TT- 


24.9 


10.3±6.9 


10.8 


2.8±1.7 


0.633±0.673 




35.6 


46.0 ±12.0 


15.1 


7.3±2.8 


0.290±0.132 


K+K--K 


24.2 


21.6±6.4 


10.9 


0.9±1.0 


0.098±0.114 


7r + 7r _ ?7(77) 


30.6 


23.7±6.9 


14.8 


0.0±1.0 a 


0.000±0.083 


7T + 7T ~ ?] ( 7T + 7T - 7T ) 


16.4 


12.7±4.8 


7.3 


1.0±1.0 


0.205±0.225 


Total 




220±22 




17.5±4.5 


0.178±0.049 
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FIG. 8: Data events (open histograms) and Monte Carlo background estimate (shaded histograms) 
of reconstructed r\ c candidate mass projection for M(ir° recoil) = 3524 ± 8 MeV/c 2 . 
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FIG. 9: Fitted 7r° recoil mass of h c candidate for M{r] c ) = 2982 ± 50 MeV/c 2 in exclusive analy- 
sis. Data events correspond to open histogram; Monte Carlo background estimate is denoted by 
shaded histogram. The signal shape is a double Gaussian, obtained from signal Monte Carlo. The 
background shape is an ARGUS function. 
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FIG. 10: Fitted photon recoil mass in data (ip(2S) — ► 77/c exclusive analysis). The signal shape is a 
double Gaussian convolved with a Breit-Wigner function. The mass resolution function is obtained 
from signal Monte Carlo. The background shape is a first-order polynomial function. The r/ c mass 
is fixed at the value 2979.7 MeV/c 2 . 
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VIII. SYSTEMATIC ERRORS 



The systematic errors on M(h c ) and B^Bh are summarized in Table IXIIII The follow- 
ing subsections describe how these errors were obtained in the individual analyses. When 
different approaches yield different results, the most conservative value is entered. 

TABLE XIII: Comparison of systematic errors in M(h c ) and B^Bh for inclusive and exclusive 
analyses. N/A: not applicable. 





M(h c ), MeV/c 2 


B 1 x B 2 


x 10 4 


Systematics in 


Inclusive Exclusive 


Inclusive Exclusive 


Number of ip(2S) 


N/A 


N/A 


0.1 


N/A 


B(^(2S) - 7%) 


N/A 


N/A 


N/A 


0.8 


Background shape 


0.3 


0.2 


0.2 


0.3 


7T° energy scale 


0.2 


0.2 


~ 


0.1 


Signal shape 


0.1 


0.1 


0.3 


0.2 


h c width 


0.1 


0.1 


0.3 


0.2 


7T° efficiency 


~ 


~ 


0.2 


0.3 


El Photon efficiency 


~ 


~ 


0.2 


0.2 


Binning, fitting range 


0.1 


0.1 


0.3 


0.2 


Modeling of h c decays 


0.1 


0.3 


0.3 


~ 


n c mass 


0.1 


0.2 


0.1 


0.1 


r/ c width 


~ 


~ 


0.2 


0.1 


rj c branching ratios 


N/A 


~ 


N/A 


0.1 


Sum in quadrature 


±0.4 


±0.5 


±0.7 


±1.0 



A. Inclusive analyses 

1. Choice of background. 

Final results in the -E(7ei) analysis were obtained using the 7r° recoil background gener- 
ated from the data. To estimate the systematic error due to choice of background, data were 
also fitted with a generic Monte Carlo background shape, yielding systematic uncertainties 
AM(h c ) ~ 0.2 MeV/c 2 , AB^B h ~ 0.2 x 10" 4 . A similar value of AB^B h was obtained in the 
M(r] c ) analysis by replacing generic Monte Carlo background by a second order polynomial 
plus an ARGUS function. However, a slightly larger value of AM(h c ) ~ 0.3 x 10~ 4 was seen 
both in data and in Monte Carlo. It is this value we quote in Table IXTTTl 

2. Photon energy calibration for 7r° energy scale. 

The standard CLEO Csl calorimeter calibration was used. To determine if the uncertainty 
in this calibration can lead to systematic error in E(jr°), the total deposited calorimeter 
energy was varied by amounts estimated by studies of radiative transitions in ip(2S) j3j| 
and 7r° — > 77 found in data. The analysis procedure, including fitting, was then repeated 
with Monte Carlo data to check for dependence on absolute calibration of CC energy. The 
small effects found may be ascribed in part to the compensating effect of the demand that 
the two photons in the low-energy 7r° have the correct effective mass. We assign an error 
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of ±0.2 MeV/c 2 in M(h c ) to the ir° energy scale on the basis of the arguments advanced in 
the subsection on the exclusive analysis. 

3. Signal shape. 

The systematic uncertainty due to uncertainty in the ir° line shape was found by varying 
the Gaussian part of the signal shape by 10% to account for a possible mis-modeling (via 
Monte Carlo) of photon energy resolution to be ~ 0.1 MeV/c 2 in M(h e ), and ~ 0.3 x 10~ 4 
in B^B h . 

4- Choice of h c resonance width. 

The systematic uncertainty due to variation of T(h c ) (0.5, 0.9, 1.5 MeV) was found to 
be ~0.1 MeV/c 2 in M(h c ), and ~ 0.3 x 10~ 4 in B^Bh- Variation of the Gaussian widths by 
±10% led to negligible changes in mass and combined branching ratio. 

5. Binning and fitting range. 

In the £"(7ei) analysis the systematic uncertainty due to fit using 1 MeV/c 2 bins, instead 
of the usual 2 MeV/c 2 bins, and changing the fitting range from 3496-3552 MeV/c 2 to 3500- 
3540 MeV/c 2 (see Table |IVJ) was found to be < 0.1 MeV/c 2 in M(h c ), and < 0.2 x 10~ 4 
in B^Bh- The M(r) c ) analysis chose 1 MeV/c 2 bins to utilize the good M(h c ) resolution 
anticipated from Monte Carlo simulations. Results were compared with those from 2 MeV/c 2 
bins and agreed with those just quoted. For the fitting range 3505-3551.2 MeV/c 2 in this 
analysis, however, B^Bh in data rose by 0.3 x 10~ 4 . This change was included as a systematic 
error associated with fitting. 

6. Modeling of h c decays. 

The signal Monte Carlo used in the E(jei) analysis took 100% of h c decaying to to ^rj c . 
An alternative signal Monte Carlo, in which 37.7% of h c were taken to decay to 7^ c and 
the rest to three gluons was generated and used to redetermine efficiency. The resulting 
B^Bh changed by ~ 0.1 x 10~ 4 . However, in the M(j] c ) analysis, larger differences were 
observed in Monte Carlo simulations when comparing B(h c — > 7?7 C ) = 37.7%, B(h c — > ggg) = 
56.8%, B(h c — > jgg) = 5.5% (nominal) and B(h c — > 7?? c ) = 100%. The nominal choice gave 
about 10% higher efficiency since events of the form ip(2S) — > 7i°h c with h c — > ggg or 
h c — > rygg sometimes pass signal selection criteria. The systematic error of 0.3 x 10~ 4 quoted 
in Table IXIIII reflects this larger value. 

7. Selected M(t] c ) range. 

In the M(t] c ) inclusive analysis, the 13 small Monte Carlo samples show that neither 
M(h c ) nor B^Bh is very sensitive to the selected M(j] c ) range in the intervals 2940-3020, 
2945-3015, 2950-3010, 2955-3005, and 2960-3000 MeV/c 2 , leading to errors of ±0.1 MeV/c 2 
in M{h c ) and ±0.1 x 10~ 4 in B^B h . 

8. Removal of "pull mass" requirement on signal 7r° . 

Instead of requiring that the signal 7r° possess the best "pull mass" within 2.5cr, all 
two-photon combinations with M(tt ) 2 within 2.5er of the correct value were considered in 
the M(i] c ) analysis. The maximum signal significance as measured by likelihood difference 
in Monte Carlo was reduced from 17.3a (Table IVjl to 16. 1<t for the nominal M(r) c ) range 
2945-3015 MeV/c 2 . Although M(h c ) obtained in the data shifted by ±0.1 MeV/c 2 from the 
nominal value, while the branching ratio shifted by ±0.9 x 10 -4 from the nominal value, 
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these shifts are within the statistical errors. No such shifts were detected in Monte Carlo 
simulations. Consequently, systematic errors were assigned to the effect of removing the 
pull mass requirement on the signal ir° of less than 0.1 MeV/c 2 in M(h c ) and 0.1 x 10~ 4 in 
B^B h . 

9. Number of neutral pions in signal region. 

Both inclusive analyses require that there be only one 7i~° candidate yielding a recoil h c 
mass within 30 MeV/c 2 of 3526 MeV/c 2 . The effect of relaxing this condition was noted. In 
all cases (independently of other selection choices), it led to Monte Carlo significances which 
decreased by 0.2-0. 3a, a decrease of M{h c ) by about 0.1 MeV/c 2 and B^Bh by 0.3 x 10~ 4 in 
data, but negligible changes in M{h c ) and B^Bh in Monte Carlo. Systematic errors in M{h c ) 
and B^pBh from this source were estimated to be less than ±0.1 MeV/c 2 and ±0.1 x 10 4 , 
respectively. 

10. Mass ranges for ip(2S) — > XJ/ip cascade suppression. 

In the M(r] c ) analysis, nominal mass ranges to suppress ir + 7T~J/if), 7r°7r°J/ijj, and j'yJ/ijj 
cascades involve recoil masses differing from M(J/ip) respectively by 8.4 MeV/c 2 (tt + 7t~), 
32 MeV/c 2 (7r°7r°), and 40 MeV/c 2 (77). These values were varied over the respective 
ranges 6.4-10.4, 22-42, and 30-50 MeV/c 2 . The maximum variations from each mode were 
then added in quadrature. Possible changes of ±0.2 MeV/c 2 in M(h c ) and ±0.2 x 10 -4 in 
B^pBh were seen in data, but negligible changes occurred in Monte Carlo simulations. These 
sources were thus estimated to lead to systematic errors of AM(h c ) < 0.1 MeV/c 2 and 
AB^B h < 0.1 x 10~ 4 . 

11. Minimum energy requirements on photons. 

In suppressing 7r°7r° J/ip cascades, a minimum energy of 50 MeV was taken for photon 
daughters in the M{j] c ) analysis. The result of reducing this energy to 40 MeV was a stronger 
suppression of both background and signal, leading to an upward shift of the mass by 0.2 
MeV/c 2 in data and no change in B^Bh in data. Changes in mass and B^Bh were negligible 
in Monte Carlo. 

12. Correction for updated M(ip(2S)) . 

The M(r] c ) analysis was based on the assumption of M{iM2S)) = 3685.96 ±0.09 MeV/c 2 , 
the world average |37| before the measurement of Ref. j2£|. With the present value of 
M(i[)(2S)) = 3686.111 ± 0.025 ± 0.009 MeV/c 2 , a correction of +0.15 MeV/c 2 thus was 
applied to the final quoted mass in that analysis. 

B. Exclusive analysis 

Because the exclusive cascade rates were measured as ratios to the radiative decays, 
systematic uncertainties related to the r\ c final state cancel. The systematic studies dealt 
with estimating the statistical significance of the h c signal, the h c mass, and the production 
branching ratio. 

In order to study the background contribution from the non--^(2S') part of the data 
(continuum data), 22 pb _1 of continuum data (beam energy ~ 1835 MeV = M(ip(2S))/2 — 
7.5 MeV) were analyzed in the same manner. The contribution of continuum data was found 
to be negligible. 

The generic Monte Carlo sample was used to see if any of the known if)(2S) decays could 
produce a fake peak which would mimic the signal. No significant peak was seen in the 
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TABLE XIV: Checks of significance, h c mass and production branching ratio (B(ip(2S) — > 7r h, 
7r°(7r/ c )) stability by varying key selection criteria (exclusive analysis). 



Selection 


Mass 


13 (cascade) / 


Significance 




(MeV/c 2 ) 


^(direct) 


M 


Default cuts 


3523.6+0.9 


0.178+0.049 


6.1 


Fit y 2 < 3 


+0.1 


0.192+0.056 


6.2 


Fit y 2 < 5 


+0.5 


0.178+0.051 


6.1 


Fit x 2 < 15 


0.0 


0.169+0.049 


5.8 


Within 30 MeV of rj c mass 


+0.7 


0.165+0.50 


5.5 


Within 40 MeV of rj c mass 


+0.2 


0.172+0.049 


5.9 


Within 60 MeV of rj c mass 


0.0 


0.172+0.049 


5.9 


Within 80 MeV of rj c mass 


-0.1 


0.188+0.052 


6.6 


Transition photon 7r° veto (2a) 


0.0 


0.168+0.051 


5.6 


Transition photon 7r° veto (4a) 


+0.2 


0.152+0.046 


5.9 


Kinematic fitted h c 


+0.4 


0.166+0.049 


5.6 


CLEOIII only 


+0.5 


0.158+0.069 


3.9 


CLEOc only 


-0.3 


0.216+0.083 


4.7 



signal region (8 bins in the 7r° recoil mass histogram, from 3516 to 3532 MeV/c 2 ) with 39 
million generic Monte Carlo events (13 times the data sample). This implies the signal seen 
in data is not due to a reflection of any known charmonium decays. 

The significance can be estimated from the background level in the signal region using 
the generic Monte Carlo or data sideband. Using events from the likelihood values of the 
fit with and without the signal contribution, we obtain s = 6.1<r; similar calculations with 
different rj c mass ranges yield s = 5.5 — 6.6a. Using events from the generic Monte Carlo 
sample, appropriately scaled so as to match event populations outside the signal region, 
we obtain an estimate of a mean background inside the signal window of 2.5 ± 0.5 events. 
Allowing for Poisson fluctuations of this number results in a probability that background 
completely accounts for the observed signal of 19 events of 1 x (s = 6.0a). The binomial 
probability that the 47 data events in Fig. and the 8 data events in the r\ c sideband, 
2600 < M(h c ) < 2860 MeV/c 2 , of Fig. El fluctuate to be greater than the 19 events in the 
signal region, 3516 < M(h c ) < 3532 MeV/c 2 , of Fig. His 2.2 x 10~ 7 , which corresponds to 
a significance of ~ 5.2a. Estimates of signal significance are summarized in Table IXIVl 

The mass of h c is estimated from a ir° recoil mass calculation. The systematic uncertainty 
associated with this estimate depends on the uncertainty of the 7r° energy scale, which is 
itself dependent on the energies of the photon daughters and their shower locations in the 
detector. Lower-energy photons and endcap photons have larger associated uncertainties. 
The fraction of endcap photons is small (<10%), so the shower- location effect on energy 
resolution was ignored. The signal ir° energy is around 160 MeV, and the corresponding tt° 
daughter photon energies vary from 30 to 130 MeV, with respective uncertainties varying 
from 1.5% to 0.2%. By changing the photon energy uniformly by ±1%, the 7r° energy in 
the signal Monte Carlo was found to shift only less than ±0.2 MeV because of the ir° mass 
constraint in the analysis algorithm which fits neutral pions. Consequently, a 0.2 MeV 
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systematic uncertainty in M(h c ) was ascribed to the 7r° energy scale. 

The t] c intrinsic width T(r] c ) has not been accurately measured. In the exclusive signal 
Monte Carlo, it is set at 27 MeV. Because the efficiency for detecting h c is estimated from 
signal Monte Carlo and a range of M(j] c ) is selected, an overestimate of V(r] c ) will result in 
an underestimated efficiency. On the other hand, it will lead to a wider signal shape for 
the f) c signal in ip(2S) — > 7?7 C and hence to an increased r) c yield. Thus the systematic error 
on the measured ratio of rates for ip(2S) — > ir°h c — > 7r°(7?7 c ) and ip{2S) — > 7?7 C is likely to 
be small because the two effects tend to cancel each other. A 2.3% systematic error was 
assigned to the ratio from the uncertainty in the rj c intrinsic width. 

The uncertainties in the rj c decay branching ratios are large; no channel is known to better 
than 25%. Changing the branching ratio of each mode 40%, once per mode, the measured 
ratio was found to shift less than 1%. Consequently, a 1% systematic error on the ratio of 
rates was ascribed to rj c decay branching ratios uncertainties. 

In the analysis of the photon recoil mass from the direct radiative decay, the r] c mass was 
fixed at 2979.7 MeV/c 2 . When this mass was floated in fitting, the value determined from the 
fit was M(j] c ) = 2970.3 ±4.1 MeV/c 2 . This result is lower than, but still consistent with, the 
CLEO inclusive photon transition study, in which the measured r\ c mass is 2976.1 ±2.3 ±3.3 
38] . Varying the fixed value of the r\ c mass in the fit of the recoil mass distribution between 
2970 and 2984 MeV/c 2 resulted in a variation of 3.6% in the yield. Half of this value, 
1.8%, was assigned to the systematic uncertainty of the combined branching ratio due to 
uncertainty in M(j] c ). 

In the decay ip(2S) — > 7i°h c — > 7c°( / yri c ), the r\ c mass selection is based on the value 
obtained by reconstructing the r\ c . When the r] c mass selection window is shifted by ± 10 
MeV/c 2 , the measured value of M(h c ) shifts by less than 0.2 MeV/c 2 . We assign 0.2 MeV/c 2 
as the h c mass systematic uncertainty due to uncertainty in M(j] c ). 

Neutral pion reconstruction efficiency has been studied in measurements of D hadronic 
branching fractions. The discrepancy between Monte Carlo and data is less than 5% [39] . 
We ascribe a 5% systematic uncertainty in the ratio of rates to 7r° efficiency uncertainty. 
This corresponds to an uncertainty in the product branching ratio of 0.27 x 10 -4 for the 
exclusive analysis and 0.18 x 10 -4 for the inclusive analysis (which finds a slightly smaller 
product branching ratio). 

In the signal Monte Carlo for the exclusive analysis, T(h c ) was set to zero, so the signal 
shape obtained from Monte Carlo essentially represented detector resolution. Varying the 
assumed value of T(h c ) up to 1.5 MeV changed the measured h c mass by less than 0.1 MeV/c 2 
and the branching ratio by 3.9%. We also studied the effects of the signal shape by changing 
detector resolution by ± 20%. The background in the exclusive study is quite small, so the 
7T° recoil mass fit range was chosen starting from 3400 MeV/c 2 . The wider background range 
helped to fit the background shape better. Varying the starting point of the fit from 3.40 to 
3480 MeV/c 2 did not change the mass and branching ratio measurement much. First- and 
second-order polynomial background shapes were used to fit the background and to study 
the systematics. The mass change was 0.2 MeV/c 2 and the rate change was 4.7%. 

The x 2 limit in kinematically constrained fits, the selection of the range for M(r) c ), and 
the veto of El transition photon candidates forming a ir° were found to be the most useful se- 
lection criteria in the exclusive study. Variation of these selection criteria within reasonable 
ranges did not change the corresponding h c mass and product branching ratios appreciably. 
The resolution in M(h c ) obtained using tt° momentum after kinematic fits was slightly better 
than that from measured E(tt°) by 2-5%, depending on modes. Because different mass reso- 



34 



lutions lead to difficulty in obtaining results and the possible gain in the mass measurement 
is small, momentum fitting was not used to obtain M(h c ). Using the kinematically fitted 
h c mass yielded values of M(h c ) and production branching ratio consistent with nominal 
results. 

IX. SUMMARY AND DISCUSSION 

A. Inclusive analyses: Summary 

Two inclusive analyses of CLEO data in search of ip{2S) — > 7t°h c — > 7r°(7?7 c ) yield an 
enhancement in the mass spectrum for recoils against 7r° attributed to the /i c (l 1 Pi) reso- 
nance of charmonium. When background is reduced by selecting a range of photon energies 
Ei'jEi) — 503 ± 35 MeV, the parameters of the resonance are found to be 

M(h c ) = [3524.4 ± 0.7 (stat) ± 0.4 (sys)] MeV/c 2 , (9) 

B^Bn = B(tp(2S) -> n°h c ) x B(h c -> 777J = [3.4 ± 1.0 (stat) ± 0.7 (sys)] x 10" 4 . (10) 

The significance of the resonance signal in this analysis, as determined by the likelihood 
method, is 3.6a. When background is reduced by selecting a range of M(rj c ) ± 35 MeV/c 2 , 
to compensate for Doppler broadening of the photon in the transition h c — > ^t\ c arising from 
the h c recoil, one finds 

M(h c ) = [3525.4 ± 0.6 (stat) ± 0.4 (sys)] MeV/c 2 , (11) 

B f B h = [3.5 ± 0.9 (stat) ± 0.7 (sys)] x 10~ 4 . (12) 
The significance of the resonance signal is 4. Oct. 

B. Exclusive analysis: Summary 

The h c produced in the reaction ip(2S) — > n°h c — ► 7T°(7^ C ) was studied by reconstructing 
■i] c in seven modes (Table UJ), leading to 17.5 ± 4.5 (stat) signal events. The significance as 
calculated from the difference in the likelihood with and without the signal contribution is 
6.1cr, and at least 5.2<r as calculated by a variety of methods. The ratio of B(ip(2S) — > 
Tr°h c — > 7r°(7?7 c )) to B{jp{2S) — > 7?7 C ) was found to be 

B(^(2S) -> 7T°h c )B(h c 7 r7 c ) , . , 

^ l L l(nQ , \ = 0-178 ± 0.049 stat ± 0.018 sys , 13 

with 

M(h c ) = [3523.6 ± 0.9 (stat) ± 0.5 (sys)] MeV/c 2 . (14) 

In CLEO III ip(2S) data, the branching ratio B(ip(2S) — ► jrj c ) was measured to be 
(3.2 ± 0.4 (stat) ± 0.6 (sys)) x 10~ 3 j3li which when combined with previous measurements 
whose average is (2.8 ± 0.6) x 1O _ ^~0, gives B(i[)(2S) -> jr) c ) = (2.96 ± 0.46) x 10" 3 . 
Combining this with Eq. ([130. one obtains a production branching ratio of 

B^Bh = [5.3 ± 1.5 (stat) ± 0.6 (internal sys) ± 0.8 (ext)] x 10" 4 , (15) 

where the last error reflects the measurement error of B(ip(2S) — > 7^ c )- The last two errors 
combine to give a total systematic error of AB^Bh = 1.0 ± 10~ 4 . 
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TABLE XV: M(h c ) and B^Bh obtained by the inclusive and exclusive analyses; combined results. 



Analysis M(h c ) (MeV/c 2 ) B^B h (units of 10" 

Inclusive E(jei) 3524.4 ± 0.7 ± 0.4 3.4 ± 1.0 ±0.7 

Inclusive M(r/ C ) 3525.4 ± 0.6 ± 0.4 3.5 ± 0.9±g;J 

Avg. Inclusive 3524.9 ± 0.7 ± 0.4 3.5 ±1.0 ±0.7 

Exclusive 3523.6 ± 0.9 ± 0.4 5.3 ±1.5 ±1.0 

Incl. + Excl. 3524.4 ± 0.6 ± 0.4 4.0 ± 0.8 ± 0.7 



C. Combination of results 

The results of the two inclusive analyses, when averaged (taking the larger systematic 
and statistical errors in each analysis), yield M(h c ) = [3524.9 ±0.7 (stat) ±0.4 (sys) MeV/c 2 
and B^Bh = [3.5 ± 1.0 (stat) ± 0.7 (sys)] x 10~ 4 . The average is taken because, as explained 
in the second-to-last paragraph of Sec. Ill, each inclusive analysis has its advantages and 
shortcomings, without a clear preference for one over the other. These results, which provide 
slightly more precise measurements of M(h c ) and B^Bh, may be combined with the exclusive 
results, based on reconstructing the 1] C in seven exclusive decay modes with much lower 
background. We have confirmed the independence of the exclusive analysis from the inclusive 
analyses by removing the exclusive signal events from our E^jei) inclusive sample. The 
results are indistinguishable from those of the original sample. We therefore combine them 
to obtain M(h c ) = [3524.4 ± 0.6 (stat) ± 0.4 (sys)] MeV/c 2 and B^B h = [4.0 ± 0.8 (stat) ± 
0.7 (sys)] x 10~ 4 , as summarized in Table IXVl 

D. Discussion 

The mass of the observed h c candidate is close to the spin-weighted average of the Xcj 
states, (3525.4 ± 0.1) MeV/c 2 . This leads to AM HF (bP) = (M(1 3 P)) - M(l 1 Pi) = [1.0 ± 
0.6 (stat) ±0.4 (sys)] MeV/c 2 , indicating little contribution of a long-range vector confining 
force or coupled-channel effects which could cause a displacement from this value. It is 
barely consistent with the (nonrelativistic) bound AM H f(1P) < |!(|. The product of the 
branching ratios for its production, B(ip{2S) — > 7i°h c ), and its decay, B{h c — > 7^ c ), is within 
the range anticipated theoretically. 
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